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摘 E 在 求解 电力 系统 经 济 分 配 问 题 时 ， 常 用 凸 二 次 函数 或 线性 函数 表示 火电 机 组 的 煤耗 成 本 。 分 段 线 
性 函数 可 以 更 精确 地 描述 机 组 煤耗 成 本 ， 但 在 求解 时 不 便于 使 用 。 为 提高 经 济 分 配 问 题 的 模型 精 
度 和 解 质量 ， 本 文 上 直接 采用 分 段 线性 凸 函数 形式 的 煤耗 成 本 模型 ， 其 次 ， 在 备用 约束 的 处 理 方 
面 ， 通 过 分 析 机 组 实际 可 提供 备用 与 机 组 实际 出 力 之 闻 非 线性 函数 关系 ， 在 不 引入 新 约束 的 前 提 
下 直接 处 理 备 用 约束 ， 最 后 ， 基 于 理论 分 析 将 模型 成 功 转化 为 线性 规划 问题 进行 求解 ， 提 出 了 两 
种 求解 方法 : 凸 组 合 系数 法 和 功率 增 莉 法 。 新 模型 和 相应 算法 在 模型 精度 和 解 质量 方面 的 性 能 均 
有 提高 ， 实 际 系统 的 算 例 测试 也 验证 了 相关 方法 的 有 效 性 。 
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1 引言 


进行 电力 系统 经 济 调度 ， 也 称 为 电力 系统 优化 调度 ， 是 电力 系统 安全 、 稳 定 和 经 济 运行 
的 要 求 。 短 期 经 济 调度 通常 以 小 时 或 若干 分 钟 为 基本 调度 时 间 单 位 确定 未 来 一 天 全 一 周 的 
机 组 组 合 和 处 于 开机 状态 机 组 的 发 电功率 ， 使 所 有 系统 约束 和 单机 组 约束 得 到 满足 ， 辐 时 
使 发 电 成 本 最 小 化 下。 短期 经 济 调度 的 数学 模型 是 一 个 大 规模 混合 整数 规划 问题 ， 且 属 
于 NP 难 (NP-Hard) 问题 ， 在 合理 时 间 内 要 得 到 全 局 最 优 解 非常 困难 。 然 而 ， 进 行经 济 调度 有 
巨大 的 潜在 经 济 效益 ， 因 此 对 其 求解 方法 的 研究 多 年 来 一 直 备 受 关注 所。 在 新 兴 的 电力 市 场 环 
境 下 ， 经 济 调度 的 意义 和 重要 性 进一步 加 强 。 例 如 ， 计 算 市 场 出 清 价 的 过 程 实质 等 价 杆 求解 一 
个 经 济 调度 问题 ;经 济 调度 目前 已 成 为 能 量 管理 系统 的 重要 子 系统 ， 对 发 电 商 而 言 ， 在 竞价 前 
或 竞价 后 进行 模拟 自 调度 或 自 调度 ( 均 对 应 于 特殊 形式 的 经 济 调度 问题 ) 是 评估 和 降低 发 电 成 
本 、 制 定 竞 价 策略 的 重要 依据 8-19。 

一 般 将 经 济 调度 问题 的 求解 划分 为 两 个 阶段 : 第 一 阶段 的 主要 任务 是 通过 求解 特定 优化 
问题 ， 获 得 各 机 纽 在 各 调度 时 段 的 启 停 状态 (离散 变量 )， 这 一 阶段 也 被 称 为 机 组 组 合 (Unit 
Commitment); 第 二 阶段 是 在 特定 机 组 组 合 下 ， 优 化 各 开机 机 组 的 发 电功率 ， 以 便 最 终 得 到 原 
调度 问题 的 一 个 近 优 解 ， 第 二 阶段 常 被 称 为 经 济 分 配 (Economic Dispatch)。 可 以 看 出 ， 机 组 
组 合 阶段 考虑 的 重点 是 离散 变量 及 其 相关 约束 ， 而 经 济 分 配 阶 段 考 虑 的 则 是 一 个 纯 连 续 优 化 
问题 。 在 某 些 算法 框架 (如 Lagrange 松 弛 ) 下 ， 第 一 阶段 和 第 二 阶段 在 算法 步骤 中 有 明确 的 分 
界 ， 另 一 些 算法 框架 (如 分 支 定 界 、 制 平面 等 一 般 的 混合 整数 规划 算法 ) 下 两 个 阶段 则 没有 明显 
Ir. 
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本 文 的 研究 重点 是 经 济 分 配 问 题 的 模型 及 求解 算法 ， 我 们 从 发 电 ) 或 发 电 商 的 角度 研究 这 
一 问题 ， 其 应 用 背景 主要 体现 在 以 下 方面 : 

首先 ， 目 前 国内 很 多 地 区 电网 中 ， 电 三 各 时 段 发 电量 指标 由 电网 预先 下 达 。 比 如 华东 电网 
在 运营 过 程 中 ， 仪 将 发 电 计划 下 达到 电厂 一 级 ， 人 允许 电 上 对 其 押 辖 机 组 进行 电量 优化 分 配 。 机 
组 组 合计 划 的 排 定 由 电网 完成 ， 电 厂 只 有 依据 电网 下 达 的 发 电量 指标 ， 对 其 开机 机 组 发 电量 进 
行经 济 分 配 的 权利 。 在 这 种 情况 下 ， 电 厂 可 以 通过 求解 经 济 分 配 问题 来 降低 煤耗 成 本 。 其 次 ， 
在 电力 市 场 环境 下 ， 当 交易 接合 完成 时 ， 机 组 组 合 、 电 价 及 各 发 电 商 分 配 的 上 网 电 景 已 定 ， 从 
而 各 发 电 商 发 电 收入 已 定 ， 但 各 电厂 仍 可 以 通过 求解 经 济 分 配 问 题 米 优化 各 机 组 出 力 以 降低 煤 
耗 从 而 增加 利润 。 最 后 ， 作 为 求解 经 济 调度 问题 两 个 阶段 中 的 一 个 阶段 ， 经 济 分 配 模型 及 算法 
性 能 对 于 提高 经 济 调度 问题 模型 精度 和 求解 效率 有 直接 影响 ， 这 也 是 本 文 研究 工作 的 主要 出 发 
点 之 一 。 

经 济 分 配 问题 一 般 是 一 个 大 规模 非 线性 连续 优化 问题 ， 目 前 求解 经 济 分 配 问 题 的 方法 主要 
包括 基于 煤耗 成 本 曲线 斜率 大 小 排序 的 优先 序 方法 、 动 态 规划 法 等 等 上 9l 。 

在 现 有 的 经 济 分 配 模型 及 算法 中 ， 有 两 个 问题 有 待 改 进 。 第 一 ， 需 要 考虑 非 线性 、 非 光 
滑 、 甚 至 非 凸 的 机 组 煤耗 成 本 模型 。 由 于 直接 考虑 非 线性 、 非 光滑 形式 的 煤耗 成 本 模型 会 给 求 
解 带 来 很 大 困难 ， 目 前 大 多 数 文献 中 考虑 经 济 调度 或 经 济 分 配 问题 时 通常 采用 凸 二 次 甚至 线性 
煤耗 成 本 模型 站， 但 在 电力 市 场 环境 下 ， 为 进一步 提高 模型 精度 ， 采 用 分 段 线性 凸 函数 形式 的 
煤耗 成 本 更 为 合适 1 ， 若 要 进一步 考虑 闪 点 效应 ， 则 必须 采用 非 凸 煤耗 成 本 模型 Q2-13] 。 

第 二 ， 现 有 模型 及 求解 方法 对 系统 备用 需求 约束 的 处 理 效率 较 低 。 文 献 14 中 ， 系 统 备 用 
约束 被 描述 为 在 任 一 时 段 所 有 开机 机 组 的 容 最 之 和 大 于 某 个 指定 值 ， 这 是 比较 粗糙 的 一 种 模 
型 。 一 种 更 好 的 表述 方式 是 在 任 一 时 段 ， 机 组 实际 可 提供 的 旋转 备用 之 和 大 于 某 个 指定 值 中， 
但 这 种 模型 在 求解 时 稍 显 麻烦 。 主 要 原因 为 : 开机 机 组 实际 可 提供 的 (10 分钟 ) 旋 转 备 用 通常 
表示 为 机 组 10 分钟 最 大 旋转 备用 、 机 组 容量 与 实际 出 力 之 差 这 二 者 中 的 最 小 值 ， 从 而 造成 机 
组 实际 可 提供 备用 与 机 组 实际 出 力 之 问 具 有 非 线性 函数 关系 ， 导 致 备用 约束 不 易 处 理 。 有 虽然 这 
种 非 线 性 函数 关系 可 以 通过 将 机 组 实际 可 提供 备用 也 作为 优化 变量 ， 同 时 引入 两 组 新 的 线性 不 
等 式 约束 来 克服 ， 但 上 述 做 法 显著 增加 了 问题 规模 ， 因 此 降低 了 求解 效率 。 

本 文 即 针对 以 上 两 个 问题 展开 研究 。 在 本 文中 我 们 直接 考虑 分 段 线 性 凸 函数 形式 的 煤耗 成 
本 模型 ， 其 次 ， 在 备用 约束 的 处 理 方面 ， 通 过 分 析 机 组 实际 可 提供 备用 与 机 组 实际 出 力 之 间 非 
线性 函数 关系 的 结构 特征 ， 在 不 引入 新 约束 的 前 提 下 直接 处 理 备用 约束 。 基 于 理论 分 析 将 经 济 
分 配 问 题 成 功 转化 为 线性 规划 问题 进行 求解 ， 提 出 了 两 种 求解 方法 : 凸 组 合 系数 法 和 功率 增 量 
法 ， 并 通过 算 例 测试 验证 了 方法 的 有 效 性 。 新 模型 和 相应 算法 在 提高 了 模型 精度 和 解 质量 的 同 
时 并 不 增加 求解 复杂 性 ， 因 此 具有 较 好 的 应 用 价值 。 

文章 结构 如 下 : 第 2 部 分 重点 介绍 本 文 提 出 的 凸 组 合 系数 法 和 功率 增 最 法 ， 第 3 部 分 是 数值 
测试 结果 分 析 ， 第 4 部 分 对 全 文 进 行 了 总 结 。 


2 经济 分 配 的 是 组合 系数 法 和 功率 增 量 法 
本 文 考虑 的 经 济 分 配 问题 的 如 下 数学 模型 


万 =min $^ C(P(t)), (1) 
Pilt) EU 
st. 3 P(t) = D(t), (2) 


icUl, 
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Y nt) > P.(, (3) 
icUl, 
P; < P(t) < Pi, icUk. (4) 


式 (1) 足 经 济 分 配 问题 的 目标 函数 。 其 中 UL 是 时 段 : 时 ， 开 机 机 组 的 集合 ， 忆 (t) 是 机 
组 i 在 时 段 t 的 发 电 功 康 ，OCi(PB(t)) 是 机 组 i 在 时 段 + 的 煤耗 成 本 ， 当 机 组 处 于 开机 状态 时 ， 
假定 共 为 分 段 线 性 是 函数 ， 如 图 1，P; 和 P; 分别 表示 机 组 i 的 最 小 和 最 大 发 电功率 ， 当 某 机 
Alio 处 于 关机 状态 时 ， 白 然 有 P.) = 0，Cio(Pio(t)) = 0。 显 然 ， 只 要 优化 开机 机 组 的 发 电 
功率 。 对 开机 机 组 i € 天 的 煤耗 成 本 曲线 ， 设 功率 区 间 上 共有 M 个 分 点 ， 分 点 (包括 两 个 端 
A A CGU (Pai, Ci) (Pas, Ci2),… , (Pam Cim)。 为 简化 记号 ， 本 文 假定 不 同 机 组 
的 煤耗 成 本 曲线 分 段 数 相同 ， 但 文中 所 有 结论 当 分 段 数 不 同 时 仍 成 立 。 式 (1) 中 的 Jg 是 经 济 分 
配 问题 的 最 优 卓 标 值 。 式 (2) 中 ，D(t) 是 电 ) 在 调度 时 段 t 内 分 配 的 总 加 发 电 功 深 ，P(t) 是 机 
组 i 的 发 电功率 ， 是 要 优化 的 决策 变量 ; 式 (3) 是 系统 旋转 备用 约束 ， 其 中 P(t) 是 时 段 t 内 系 
统 总 备用 需求 ，ri(t) 是 机 组 i 在 时 段 t 提 供 的 备用 ， 计 算 公式 为 


ri(t) = min { 7i, P; = P(t) 上 (5) 


这 里 元 是 机 组 i 能 提供 的 圾 大 备用 ， 是 一 个 常数 。 


P0 
o 2 4 6 8 10 
P, ,=Pmin Pio Pim Pi M-1 P, 7" Pmàx 


图 1: 开机 机 组 i 的 煤耗 成 本 曲线 Ci(P(t))， 这 里 Puis = P, 


需要 说 明 的 是 ， 模 型 (1)-(4) 中 虽 未 明确 考虑 爬升 约束 ， 但 爬升 约束 事实 上 是 隐 含 在 其 
中 的 ， 对 此 问题 解释 如 下 : 经 济 分 配 一 般 是 迎 时 段 进 行 的 ， 如 果 机 组 i 存在 息 升 约束 ， 不 
Jj; WEAR AE [P — Pt- 1) € Ai。 式 中 AA; 为 机 组 i 在 相 邻 两 个 时 段 发 电 功 梁 域 
大 允许 变化 量 ， 即 最 大 允许 爬升 量 。 在 求解 第 t+ 时 段 的 经 济 分 配 问 题 时 ， 第 上 -工时 段 的 经 
济 分 配 已 完成 ， 从 而 五 全 -了 J 取 值 已 定 ， 因 此 上 述 肘 升 约束 事 实 土 转化 为 关于 五 父 的 竹 约 
R P((t—1)— A; x P(t) € P(t 一 1) 十 Ai。 人 在 求解 时 ， 上 述 界 约束 和 机 组 i 的 最 小 及 最 大 发 电 
功 涩 界 约 束 一 起 被 用 村 确定 P.) 的 上 下 界 ， 即 界 约 束 的 准确 形式 为 


max {Pi, P(t—1)—Ai} € P(t) € min {P(t—1)+Ai, Pi), (6) 


zd 


在 实际 求解 时 我 们 采用 的 是 约束 (6)， 模 型 表述 中 采用 约束 (4) 只 是 为 了 形式 上 的 简洁 ， 显 
然 (4) 和 (6) 有 完全 类 似 的 结构 。 此 外 ， 可 能 出 现 的 另 一 种 情况 是 在 机 组 组 合 阶段 已 指定 在 未 来 


68 I 2 数 学 学 R 第 27 卷 


的 某 个 时 段 ， 机 组 发 电功率 被 限制 在 某 个 水 平 或 某 个 范围 内 ， 如 某 些 机 组 在 开机 区 间 的 最 后 一 
个 时 段 要 求 发 电功率 为 最 小 ， 此 时 可 按 和 上 土 述 做 法 类 似 的 方法 ， 反 方向 逐 时 段 递 推 至 时 段 电 
获得 在 时 段 + 内 人 允许 的 出 力 范 围 ， 以 上 在 考虑 爬升 约束 时 确定 PE) E FARIEY 
前 顾 后 ”。 

式 (5) 表明， 机 组 i 在 时 段 + 能 提供 的 旋转 备用 ri(t) 可 以 表示 为 机 组 发 电功率 忆 () 的 一 元 函 
数 ， 而 且 是 仅 有 一 个 转折 点 的 分 段 线 性 凹 函 数 ， 其 图 像 如 图 2 所 示 。 由 于 机 组 煤耗 成 本 曲线 也 
是 P(t) 的 一 元 分 段 线性 函数 ， 所 以 不 妨 假设 图 2 中 唯一 的 一 个 转折 点 横 坐 标 也 是 CLCPL 0) 的 
某 个 转折 点 的 横 坐 标 (否则 增加 一 个 分 点 即 可 )， 并 设 其 转折 点 编号 为 ri。 然后 进一步 将 旋转 
备用 曲线 在 形式 上 表示 为 有 《个 分 点 的 分 段 线性 函数 ， 且 所 有 分 点 的 横 坐 标 与 机 组 煤耗 成 本 
曲线 的 分 点 横 坐 标 全 部 重合 ， 因 此 旋转 备用 曲线 的 所 有 分 点 可 表示 为 


(已 ri) (Pia ri2) (Pun rM) (7) 
且 有 以 下 关系 式 成 立 
Tin = Fi, 车 1l<n< mi, 
E (8) 
Tin = Pi — Pin, d$ mi Ent M. 


以 上 处 理 将 对 理论 分 析 和 数值 计算 带 来 很 多 方便 。 

既然 关切 可 以 表示 为 机 组 发 电功率 瑟 雪 的 一 元 函数 ， 则 任意 给 出 一 组 瑚 (办 即 可 得 到 一 
组 ri(t)， 然 后 将 PO) RI ri (0) RA (2)-(4) 式 即 可 知 它们 是 否 构 成 问题 (1)-(4) 的 一 组 可 行 解 ， 
换 旬 话说， 一 组 解 是 否 可 行 完全 由 P(t) 的 取 值 确定 ， 因 此 在 后 文中 讨论 解 的 可 行 性 时 全 部 围 
R P(t) 的 取 值 进行 分 析 。 基 于 上 述 理 由 ， 在 本 文中 未 将 ri(t) 作为 经 济 分 配 问题 中 的 独立 变 


里 


图 2: 开机 机 组 i 的 旋转 备用 曲线 ri(t)， 这 里 Puis = P, 


本 文 提出 的 第 一 种 经 济 分 配方 法 一 凸 组 合 系数 法 的 主要 思想 是 充分 利用 煤耗 成 本 曲线 的 凸 
性 和 机 组 备用 曲线 的 四 性 ， 建 立 一 种 线性 规划 模型 ， 进 而 获得 经 济 调度 问题 的 最 优 解 ， 下 面 先 
给 出 一 个 引 理 。 

引 理 1 对 经 济 分 配 问 题 (1)-(4) 的 任 一 组 可 行 解 P(t) (i € UI). 一 定 存在 一 组 实数 aim (i € 
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UJ, 1 < m < M) 使 以 下 表达 式 成 立 


M 
P(t) = X din Pim, (9) 
m=i 
M 
Dum=1, 0<aim <1, (10) 
m=1 
M M 
5 Y oonPin cm D 5 S diari S P. (11) 
i€U P, m=1 i€U I, m=1 


有 反之， 车 一 组 实数 oim (i e UR, 1 <m < M) 使 (10) 和 (11) 式 成 立 ， 则 (9) 式 中 定义 
的 五 的 给 出 问题 (1)-(4) 的 一 组 可 行 解 。 

证 明 首先 ， 如 果 瑟 区 对 应 于 问题 (1)-(4) 的 一 组 可 行 解 ， 则 Bi(t) 一 定位 于 煤耗 成 本 曲线 
众多 分 点 中 某 两 个 相 邻 分 点 的 横 坐 标 之 间 ， 即 存在 m (1 < ni < M 一 1)， 使 得 下 式 成 六 


Pin; € Pilt) € Pina. (12) 


由 于 分 段 线性 的 旋转 备用 曲线 具有 和 煤耗 成 本 曲线 完全 相同 的 分 点 横 坐 标 ， 再 根据 线性 函数 的 
性 质 ， 设 解 Pi) 对 应 的 备用 为 ri(t)， 则 存在 实数 入 ; € [0,1] 使 得 


P(t) = Xi Pan; + (1 — A) Panini; (13) 
ri(t) = Ai rin; 十 (1 一 Xi)7rini+l. (14) 
按 下 述 方式 定义 aam 
0， 如 果 m<ni E m»ni-1, 
Qim = 4 A. 如 果 m= ni, (15) 


1 一 X， 如 果 m=n;+1. 


则 由 (13) 和 (15) 可 得 (9) 式 ， 由 (15) 可 得 (10) 式 ，(2)、(3)、(13)-(15) 结合 起 来 即 得 (11) XX. 
反之 ， 若 一 组 实数 aim (i EeE UI,1<m<M) 使 (10) 和 (11) 式 成 立 ， 则 对 固定 的 i，(10) 式 
表明 aim( < m < M) 恰 好 构成 一 组 凸 组 合 系数 ， 因 此 若 根据 aizm 由 (9) 式 定义 出 一 
组 P(t)， 则 P(t) 是 机 组 i 的 煤耗 成 本 曲线 及 备用 曲线 上 各 分 点 横 坐 标的 一 个 凸 组 合 ， 
因此 已 (#) 的 值 必 介 于 各 分 点 横 坐 标的 最 大 值 与 最 小 值 之 间 ， 从 而 约束 (4) 自然 满足 ， 再 
由 (9) 式 及 (11) 式 中 的 第 一 个 式 子 即 得 (2) 式 成 立 ， 因 此 仅 需 再 证 明 (3) 式 成 立 。 事 实 上 ， 如 
果 瑟 人 由 (9) 式 定义 ， 则 因为 mm 人 是 五 的 的 上 四 函数 ， 所 以 根据 上 四 函数 的 性 质 : 若干 点 的 凸 组 
合 处 的 函数 值 大 于 等 于 这 些 点 处 函数 值 的 相应 凸 组 合 ， 这 些 分 点 处 的 备用 值 为 mim， 从 而 有 


M 
nz autos (16) 
icUI, m=1 
(16) 式 和 (11) 式 中 的 第 一 个 式 子 结合 起 来 即 证 明了 (3) X. 证 毕 


引 理 1 表明 经 济 分 配 问题 的 任 一 组 可 行 解严 (人 (GE VE) 可 以 用 一 组 恰当 的 凸 组 合 系 
JE ois (i e UL,1<m<M) 米 表示 ， 同 时 引 理 1 的 证 明 过 程 也 上 暗示 这 种 表示 方式 可 能 不 唯 
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一 。 比 如 若 某 个 机 组 ;的 煤耗 成 本 曲线 分 点 横 坐 标 为 : 10、20、30， 则 同一 个 数值 15 可 以 有 以 
下 多 种 凸 组 合 表示 方式 
15=0.5x10+0.5x20+0x30=0.75x10+0x20+0.25x30， (17) 


这 种 不 唯一 的 表示 方式 给 成 本 计算 带 来 了 困难 ， 因 为 在 用 凸 组 合 方式 计算 分 段 线性 的 发 电 成 本 
时 实际 上 只 有 机 组 发 电功率 所 处 的 小 区 间 的 两 个 端点 处 的 凸 组 合 系 数 可 能 不 为 0， 其 余 分 点 的 
凸 组 合 系数 必须 为 0， 以 上 困难 通过 下 述 定 理 得 到 解决 。 

定理 1 ( 凸 组 合 系数 法 ) 考虑 如 下 线性 规划 问题 


M 
(LP1) Jf = min » DD Qim Cim , (18) 
tiom CUI, m-i 
变量 所 受 约束 为 (10)-(11) 式 ， 民 为 最 优 目标 函数 值 ， 以 om € UI, 1€ m € M) 
示 (LP1) 的 任 一 组 最 优 解 ， 则 按 下 述 (19) 式 得 到 的 Pr(t) (i € UJ) 是 问题 (1)-(4) 的 一 组 最 
优 解 。 
PEO = > im Pn. (19) 


证 明 首先 证 明 Jf — Jd. 

对 于 (LP1) 的 任 一 组 可 行 解 aim ( € UB, 1 <m<M)， 由 引 理 1 可 知 车 按 (9) 式 计算 出 一 
2H P(t). BU P,(t) teih EE (1)-(4) 的 一 组 可 行 解 。 注 意 到 煤耗 成 本 曲线 是 下 凸 的 ， 据 凸 函数 的 
性 质 : 若干 点 的 凸 组 合 处 的 函数 值 小 于 等 于 这 些 点 处 函数 值 的 相应 凸 组 合 ， 从 而 有 


M 
Ci(P(t)) < »» Oi,m Cim - (20) 
m=1 
于 是 
M 
2.9(B(00)S Y, J, aw (21) 
EUH i€UI, m=1 


因为 P;(t) (i € U I) 4 i8) 8 (1)-(4) 的 一 组 可 行 解 ， 所 以 (21) 式 的 左边 是 问题 (1)-(4) 的 可 行 值 ， 
因而 


J< $ (P0). (22) 
icUl, 
Hi (21) & (22). # 
M 
AE YO 6m Cn: (23) 
i€UlI, m=1 


对 比 (23) 式 的 右边 与 问题 (LP1)， 而 Qim (i € Ule, 1 < m < M) 是 问题 (LP1) 的 任 一 组 可 行 
解 ， 于 是 
Jas. (24) 
KTR, W P(t) (i e UL) 是 问题 (1)-(4) 的 一 组 最 优 解 ， 由 引 理 1， 一 定 存在 问题 (LP1) 的 
一 组 可 行 解 Qim (i € UL, 1 < m< M)( 按 (15) 式 中 的 取 法 )， 必 然 有 


C(BR.(0) = 》 aim Cun. (25) 
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Jg — 5 C;(P.(t)) = 5 Tain im (26) 


iEU r icUI, m=1 
注意 (26) 式 中 第 一 个 等 式 成 立 是 因为 现在 忆 (t) (i € UI) A eL f (1)-(4) 的 一 组 坡 优 解 。 
而 aim (iE UJ, 1<m<M) 是 问题 (LP1) 可行 解 ， 于 是 


EST (27) 


(24) 和 (27) 结合 起 来 即 证 明了 Jr = Jo. 

其 次 ， 证 明 (19) 式 给 出 问题 (1)-(4) 的 一 组 最 优 解 。 

根据 给 定 条 件 ， 设 ar ， (ie UL, 1I<m<M) 是 (LP1) 的 任 一 纽 最 优 解 ， 则 因为 (19) 式 实 
际 与 (9) 式 结构 相同 ， MO BEL, Fà (19) 式 得 到 的 P* (t) (i € UI) t& Ae fe] E (1)-(4) 的 一 组 可 
行 解 ， 基 于 和 本 定理 证 明 中 (21) 式 以 前 的 部 分 完全 相同 的 推 埋 可 得 


M 
TAR Y, Y anO (28) 


i€Ut, i€EUT: m=1 


注意 (28) 中 的 Px(t) 是 问题 (1)-(4) 的 一 组 可 行 解 ， 而 at 是 (LP1) 的 一 组 最 优 解 ， 因 此 根据 两 
个 问题 最 优 目 标 值 的 相关 符号 ，(28) 式 可 扩充 为 下 式 


M 
As c0 s9) 60-4. (29) 

iEUT: icUI, m=1 

因为 我 们 已 证 明了 Jr = 大， 所 以 (29) 中 的 所 有 不 等 号 全 部 变 为 等 号 ， 从 而 
= M eR). (30) 
EU 

(30) 式 表明 问题 (1)-(4) 的 可 行 解 P (£) (i € UE) 对 应 的 目标 函数 值 等 于 最 优 值 ， 所 以 它 必 然 是 
一 组 最 优 解 。 证 毕 


下 询 从 另 一 种 角度 将 经 济 分 配 问题 转化 为 一 个 规模 比 (LP1) 更 小 的 线性 规划 问题 ， 为 此 先 
对 有 关 的 两 个 符号 加 以 解释 


. Rx Cin4i un Cin 
Pin B. 一 了 


, n—1,2,..,M-—1, (31) 
indi — Pin 


ERIALA i MIRREN Ha x LI n BRRRER, 常量 。 对 于 经 济 分 配 问 题 的 任 一 个 可 行 
解 ， 机 组 i 的 发 电功率 P(t) dcc TS BL, WR 


Pim < P(t) < Pi,m+1: (32). 
XpEASm(1£mtzM -1 R, 则 按 如 下 方式 定义 一 组 变量 APnGeUn, 1 Em € M — 1) 
APn=Pntii-Pn, d )izntzm, 
AP; n = O0, di m«nzM-1 
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称 APin 为 功率 增 量 。 容 易 验 证 ， 如 果 机 组 i 的 发 电功率 P(t) 不 与 煤耗 成 本 曲线 上 任 一 个 分 
点 的 横 坐 标 重合 ， 则 (32) 中 的 m 是 唯一 的 ， 因 此 (33) 式 唯一 地 确定 了 功率 增 最 ， 如果 P;(t) 5 
某 一 分 点 的 横 坐 标 重合 ， 则 (32) 中 的 m 不 唯一 ， 但 对 (33) 式 定义 的 功率 增 景 取 值 无 影响 ， 
即 (33) 式 定义 的 功率 增 景 取 值 始 终 是 唯一 的 。 
将 (33) 式 中 所 有 等 式 全 部 相 加 ， 立 即 得 到 
M-1 M-1 
P(t) = Pi 》 AB m P; Ñ AB. (34) 
n=1 n-l 
(33) 和 (34) 结合 起 来 表明 由 机 组 发 电功率 可 以 唯一 确定 功率 增 量 ， 反 过 来 由 功率 增 量 也 可 以 唯 
一 确定 机 组 发 电功率 。 因 此 ， 在 本 文 提出 的 将 经 济 调度 问题 转换 为 线性 规划 的 第 二 种 方法 中 ， 
直接 考虑 的 是 如 何 得 到 功率 增 量 ， 我 们 有 下 面 的 定理 。 
定理 2 (功率 增 量 法 ) 考虑 如 下 线性 规划 问题 


M-1 
(P2) B= min $^ Y, pin APim (35) 
i€EUT: n=1 
st. 0< APi;,n < AP;,n+1 us Pin, tE Url, n = 1,2,- is ,M = l, (36) 
M-1 
$3 XAaBa-D()- Y, B. Gn) 
icUI, n=1 iEUT: 
M-1 
Y DAPn < 》 n-R() (38) 
i€U I, n-mi iEUT: 


其 中 (38) 式 中 的 mi 在 (7)、(8) 式 及 其 前 面 一 段 文 字 中 已 定 义 过 了 。 问 题 (LP2) 与 原 经 济 分 配 
问题 (1)-(4) 有 如 下 关系 。 

(a) 3$ P(t) (i € UE) 是 问题 (1)-(4) 的 一 组 最 优 解 ， 则 按 此 忆 (t) 根 据 (33) 式 得 到 的 功率 增 
量 是 问题 (LP2) 的 一 组 最 优 解 。 

(b FAPinli EUl, Il<mn<M-1T) 是 问题 (LP2) 的 一 组 最 优 解 ， 则 按 此 解 根据 (34) 式 
得 到 的 P; (t) (i € U I5) 是 经 济 分 配 问题 (1)-(4) 的 一 组 最 优 解 。 

证 明 证 明 分 三 步 。 

第 一 步 ， 首 先 证 明 若 P(t) (i e UL) 是 问题 (1)-(4) 的 一 组 最 优 解 ， 则 按 此 P(t) 根据 (33) 式 
得 到 的 功率 增 量 是 问题 (LP2) 的 一 组 可 行 解 。 

如 果 已 人 (人 € UH) 是 问题 (1)-(4) 的 一 组 最 优 解 ， 则 对 比 (32)、(33) 与 (36) 式 ， 可 知 按 
此 P(t) 根据 (33) 式 得 到 的 功率 增 量 APin (i € UL 1 < n € M1) 必然 满 足 约 束 (36)， 如 前 
所 述 ， 由 (33) 可 得 (34)， 而 约束 (37) 是 (2) 式 与 (34) 式 的 直接 推论 ， 如 果 再 证 明了 AP, 还 满 
足 (38)， 则 上 述 分 析 表 明 根 据 (33) 式 得 到 的 功率 增 量 AP, n 是 (LP2) 的 一 组 可 行 解 。 事 实 上 ， 
根据 (5) X m, 的 定义 可 得 


(39) 


Ti = 


Ti, 若 P(t) < Pima 
Pi- P(t), E PPB(t) > Pimi, 


Pj = Fi + Pimi- (40) 
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再 根据 (33) 式 可 知 
M 0, 若 P;(t) < Pimi» 
Tm P(t) Pim; 若 P(t) 2 P; mi- 
Ki (40). (41) 代入 (39) 立即 可 得 
M-—1 
ri(t) = Ti — b» APin. (42) 


n—mi 


(3) RAI (42) 式 结合 起 来 即 证 明了 (38) 式 。 现 在 我 们 已 证 明了 AP 是 (LP2) 的 一 组 可 行 解 ， 
如 果 记 该 可 行 解 对 应 的 目标 值 为 及， 则 根据 有 关 记 号 、 分 段 线性 函数 的 性 质 以 及 P(t) (1 < i < 
T) 是 问题 (1)-(4) 的 最 优 解 这 个 假设 立即 得 到 


M-1 
$3 APRO) = > Ci(P)+ 3, Y, sis AB 


Jg = 
iEUL, iEUT: ieUT: n=1 
= M G(P)tA2 》 GP) + J3. (43) 
i€UI, icUl 


58 b. KEHE AP, nl € UL,1€nExM 1) (LP2) 的 一 组 最 优 解 ， 则 按 此 
解 根据 (34) 式 得 到 的 P(t) (i e UI, 1 < i < 了 是 经 济 分 配 问题 (1)-(4) 的 一 组 可 行 解 。 

我 们 先 对 解 APin Ge Uk,1<n<M-1) 进 行 必要 的 技术 处 理 。 考 虑 机 组 i 在 煤耗 成 本 曲 
线 各 个 分 段 的 功率 增 量 AP, APi,2,… , AAP,m-1， 如 果 存 在 某 两 个 段 n1, n 使 得 下 述 关 系 
成 立 

^i«ns, H APn — Pntil, AB. (44) 


即 编 号 较 小 的 段 的 功率 增 量 未 达到 最 大 ， 但 编号 较 大 的 段 的 功率 增 量 非 0， 则 我 们 可 以 适当 
增加 入 Pi 同时 减少 人 Pn, (增加 与 减少 的 量 相 同 ) EE AP, S, 增加 人 至 其 上 限 或 APin, 减少 
至 0 以 确保 (36) 式 始终 满足 ， 以 上 处 理 显然 不 会 改变 功率 增 量 AP, APi2z,… ,APm-i 之 
和 ， 所 以 (37) 式 始终 成 立 ， 对 于 (38) 式 ， 以 上 处 理 只 会 使 得 不 等 号 左边 更 小 ， 从 而 确保 经 上 述 
调整 后 仍 得 到 (LP2) 的 可 行 解 ， 但 由 于 煤耗 成 本 曲线 是 下 凸 的 ， 所 以 各 段 斜 率 pin Bl m RK 
而 增 大 ， 所 以 对 比 里 标 函 数 表达 式 (35) 可 知 经 上 述 调整 后 目标 函数 值 不 会 增加 ， 又 因为 已 假 
设 APn(ieEUL,1<n< M1) 是 问题 (LP2) 的 一 组 最 优 解 ， 所 以 经 上 述 调 整 后 得 到 的 始终 
是 (LP2) 的 最 优 解 ， 现 在 将 上 述 调整 过 程 一 直 进行 全 不 能 再 调整 为 止 ， 记 最 后 的 解 为 


AB, ieUlh, 1€nxM-1. (45) 


根据 以 上 分 析 我 们 可 以 得 知 功率 增 量 A (i EU, 1n E M - 1) RGULFTEUR. 

1) APBn(ieEUL,1<n< M1) 是 (LP2) 的 最 优 解 。 

2) 由 AB 代入 (34) 式 得 到 的 P(t) 与 由 APin 代 入 (34) 式 得 到 的 忆 (t) 相同 。 

3) 如 果 某 个 AP; > 0， 则 当 7 < n 时 AP,; = P; j+1 — Ho 如 果 某 个 AP, « Pin4A = 
Pin: Jj 24 j > nit AB; =Q. 

MEZE H AP, n 代入 (34) REEI P(t) (Eh AP; n 代入 (34) 式 得 到 的 P(t), (36) 和 (34) 
结合 起 来 保证 了 (4) 式 成 立 ，(34) RI (37) 结合 起 来 保证 了 (2) 式 成 立 ， 利 用 得 到 的 P(t) 可 进 一 
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步 计算 出 7i(t)， 仿 照 定理 2 证 明 中 第 一 部 分 (42) 式 以 前 的 分 析 推 理 可 得 


M-1 " 
rt) =T- 9, AP. (46) 
再 根据 AP, n 与 APin 的 关系 可 得 
M-i . M-1 
ri(t) =ñ; 一 > AP; n > 7i 一 »» APia. (47) 


(47) XJ (38) 式 结合 起 来 即 得 到 (3) 式 ， 从 而 书信 (L<i< 站 是 经 济 分 配 问题 (1)-(4) 的 一 组 
可 行 解 ， 记 该 可 行 解 对 应 的 目标 函数 值 为 nn， 则 根据 以 上 分 析 及 有 关 记 号 有 下 述 结论 成 立 


Jg € Jo 5 C; (P.(t)) 


icUl, 


M-1 
35 60) 》 Moss AB. = Y OQ. (48) 


i€EUT: i€EUI: n-l icUI 


第 三 步 ， 最 后 证 明 结 论 (a). (b) 成 立 。 

对 比 (43) RH (48) 式 可 知 两 式 中 所 有 不 等 式 实际 上 成 为 等 式 ， 从 而 大 = Jz, Jo = Jg. BH 
结论 (a) (b) 成 立 。 证 毕 

注 1 仔细 观察 经 济 分 配 问 题 的 数学 模型 (1)-(4) 式 可 以 发 现 ， 目 标 函 数 (1) 及 约束 (3) 中 
均 含 有 非 线性 因素 ， 因 此 直接 求解 会 比较 麻烦 ， 定 理 1 和 定理 2 的 意义 在 于 充分 利用 经 济 
分 配 问题 的 结构 特性 将 其 转化 为 简单 的 线性 规划 问题 求解 ， 因 而 具有 重要 意义 。 归 外 ， 在 问 
BL(LP2)'h, 358 (M—1)|U L| 25&, 2(M-1)-UD|-24 2498. LPE 2(M —1)-U I| 
界 约束 ， 因 此 变量 和 约束 个 数 均 少 于 问题 LPi， 但 需要 根据 (31) 式 计 算 (M — 1) ULAR 
率 值 。 

注 2 ”本文 讨 论 的 仅 是 单 时 段 的 最 优 经 济 分 配 问 题 ， 当 考虑 多 时 段 经 济 分 配 时 ， 可 以 将 本 文 
提出 的 方法 按时 段 兴 代 应 用 杆 各 时 段 并 最 终 得 到 多 时 段 问 题 的 解 。 男 一 种 可 能 途径 是 直接 整体 
求解 多 时 段 经 济 分 配 问 题 ， 此 时 间 题 规模 将 急剧 增 大 ， 特 别 是 当 考 虑 网 络 传输 安全 约束 时 问题 
复杂 性 进一步 增加 。 目 前 ， 我 们 已 尝试 将 本 文 提出 方法 扩展 全 考虑 安全 约束 的 多 时 段 经 济 分 配 
问题 整体 求解 算法 设计 中 ， 有 关 研 究 正在 进行 。 

注 3 本 文 用 分 段 线性 函数 直接 近似 二 次 成 本 函数 的 做 法 是 有 意义 的 ，2009 年 翟 桥 杜 等 在 
文献 [16] 中 给 出 了 分 段 线性 成 本 函数 代替 -次 成 本 函数 时 ， 其 误 芜 依赖 于 等 分 区 问 的 长 度 。 


li 


3 测试 算 例 


我 们 采用 IEEE 标 准 测试 数据 中 测试 我 们 的 算法 ， 这 是 一 个 26 台 机 纽 的 电力 系统 ， 机 
组 :的 煤耗 成 本 曲线 用 分 为 20 段 的 分 段 线性 凸 函数 表示 ， 其 对 应 的 凸 二 次 曲线 逼近 形式 为 


Ci(Pi(t)) = a.P?(t) + biP;(t) + ci, 


机 组 出 力 上 下 限 、 煤 耗 曲线 数据 系统 及 各 小 时 功率 需求 表 参 见 文献 15]。 系 统 的 旋转 条 
H P.(t) 79 RB ok D(t) B 576. 
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计算 中 ， 我 们 先 采 用 Lagrange 松 弛 法 求解 机 组 组 合 问题 ， 得 到 一 个 近 优 的 机 组 组 合 后 ， 再 
采用 本 文 提出 的 方法 求解 经 济 分 配 问题 。 由 Lagrange 松 弛 法 得 到 的 机 组 组 合 结果 在 表 1, 中 给 
出 ， 共 中 “13” 表 杰 机 组 在 对 应 时 段 处 于 开机 状态 ，“0” 表 杰 处 于 关机 状态 。 


表 1: X Lagrange 松弛 方法 获得 的 可 行 机 组 组 合 


du 每 个 调度 时 段 的 机 组 状态 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 
1 |1 0 0 0 0 0 O 1 11 00000001 10 1 0 0 0 
2 |1 0 0 0.0 00 1 1 ı1 0 0 0 0 0 0 O 1 1.0 10 0 0 
3 |1.0 0 00 00 1 1 1 0 0 0 0 0 0 O 1 1.0 10 0 0 
4 |1 0 0 0 0 0 O 1 1 1 0 0 0 0 0 0 O 1 1 0 10 0 0 
5 |1 0 00000 11 ı1 0 0 0 0 0O O 0 1 1 0 100 0 
6 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
7 lo o o o o 0o o o 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
8 |o o o o o o o o 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
9 10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
10 |1 1 1 1 1 1 1 1 1 1 1 10 1 1 10 1 1 1 1 1 1 1 1. 1 
11 1 1 0 0 1 1) 10 140 14 10 10 1) 1) 1) 1 1) 1 1 1 1 14 10 1. 1 
12 |1 1 0 0 0 0 1) 1 1 10 1 1 1) 1 1 1 1 1 13 1 1 1 1 1 
13 |1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
14 |1 1 1 1 0 0 0 0 1 1 1 1 1) 1 1 1 1 1 1 1] 1 1 1. 1 
15 |1 1 1 10 0 0. 0 0 1) 10 1) 109 1 1 1 1 1) 12 1 1.) 1 1. 1] 
1610 0.00 0.0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
ipn ub. b C) dod do dd rl S oL.4 4.1 1 1-c1.l1 l1 ro i d 
18 |1 1 1 1 10 10 1 1 1 1 1 10 1 1 10 1 1 1 1 1 1.) 1 1 1 
19 |1 1 10 1 ) 1 .) 1) 1 1) 1 1 10 1 1 1 10 1 1 1) 1 1.) 1 1 
20 |1 1 1 1 1 1 1) 1) 1 1) 1 10 1 1 1 1 1 1 1 1.) 1 1.) 1. 1 
2 |o 0 0o 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
2 |o 0 0 0 0 0 0 0 O0 1 1 1 1 1 1 1 1 1 1 1 1 1) 1 0 
238 jo 0 0 00 0000000 00 0 0 1 1 1 1 1 1 0 0 
24 |1 1 1 10 10 1 1 1 1 1 1.1 1 10 1 1 1 1 1 1 1 1 1 1 1 

1 1 1 10 1 1 1 1 1 1 1 1 11 1 1 1 1 1 1 1 1 1 1l 1 

P xo dod dob do dL d Y i1 Ilo 1 Y 1 1 d YT 4 1 


其 次 ， 在 每 一 个 调度 时 段 ， 依 照 文献 中 ， 将 机 组 分 为 发 电功率 可 调 机 组 和 不 可 调 机 组 ， 然 
后 ， 采 用 本 文 提出 的 ， 得 到 每 个 调度 时 段 内 每 台 可 调 机 组 的 煤耗 成 本 曲线 对 应 的 组 合 系数 (或 
功率 增 量 )， 按 照 公 式 (9) (或 者 公式 (34)) 和 公式 (5) (或 公式 (42)) 得 到 可 调 机 组 在 指定 时 段 的 
发 出 功率 和 和 备用， 而 其 它 机 组 的 发 电功率 和 备用 保持 为 对 偶 解 的 数值 不 变 。 最 后 ， 即 得 到 经 济 
分 配 问题 (1)-(4) 的 最 优 解 ， 将 经 济 分 配 的 解 与 机 组 组 合 结果 结合 起 来 ， 即 得 到 了 整个 经 济 调度 
问题 的 最 终 解 。 

在 第 一 时 段 ， 可 调 机 组 是 1-6，10-15，17-20，24-26， 共 有 19 台 功 率 可 调 机 组 。 由 于 篇 幅 
原因 ， 仪 列 出 第 一 小 时 各 台 机 组 的 发 电功率 PA) 和 备用 值 7;(1)， 见 表 2。 计 算得 到 的 24 小 时 
对 偶 成 本 739329.88 美 元 ， 总 发 电费 用 为 743128.5 美 元 。 对 偶 间 际 0.51%%， 由 于 经 济 分 配 的 结 
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果 对 最 终 解 质量 也 有 重要 影响 ， 如 此 小 的 对 偶 间 隙 恰好 验证 了 本 文采 用 的 经 济 分 配 模型 及 相应 
的 求解 算法 是 有 效 的 。 


表 2: 第 一 小 时 各 台 机 组 的 发 电功率 和 备用 


ra | Roa | «cow [a | non [ «ron 
princ qc Cur EIE 
£o a 39 [5s] 9 v 
ipu cs Spec]: 1L 9 
C] a | o3 que] o | o 
DE | o9 qom | oss 
[04 pow [q»9[ s 
r| o | o fo x [| 
| o Ta。 | , 
8| 9 | |] * ] 3 
e| mw | o fa] o> |` 
u | e | o |a] w i 
u| mw | » |s) 4 | ^? 
s w | o | w i 


对 经 济 调度 和 经 济 分 配 问题 ， 分 段 线性 函数 形式 的 煤耗 成 本 模型 更 为 精确 但 通常 会 给 求解 
带 来 不 便 。 如 果 采 用 分 段 线性 凸 函数 形式 的 煤耗 成 本 模型 ， 则 结合 备用 与 机 组 出 力 间 的 四 函数 
关系 可 将 经 济 分 配 问题 转化 为 线性 规划 问题 求解 ， 从 而 可 在 模型 精度 和 模型 易 解 性 方面 同时 满 
足 需求 ， 实 际 系统 的 算 例 测 试验 证 了 相关 方法 的 有 效 性 。 
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New Methods for Economic Dispatch in Power Systems Based on 
Linear Programs 


GUO San-gang, ZHANG Lin, CAO Jii, LI Xiao-kang 


(Department of Mathematics, Shaanxi University of Technology, Hanzhong 723001) 


Abstract: The fuel cost of thermal units is often approximated by linear or convex quadratic functions 
when solving the economic dispatch (ED) problems. Piecewise linear functions provide more precise 
formulation of fuel costs whereas make the ED problems much harder to solve. A convex piecewise 
linear fuel cost formulation is adopted in this paper to improve the accuracy of the model. The reserve 
constraints are re-formulated as constraints on generation levels of units, based on a detailed analysis 
of the structure of the constraints. Then, the ED problem with a convex piecewise linear fuel cost 
formulation is transformed into two kinds of linear programming (LP) problems: an LP formulation 
based on convex combinatorial coefficients and an LP formulation based on power output increments. 
The new methods are finally implemented, and numerical results of real systems show that our methods 
are efficient and effective. 

Keywords: power system; economic scheduling; economic dispatch; linear programming 
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